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The Fat Link Irrelevant Clover (FLIC) fermion action is a variant of the 0(a)-improved Wilson 
action where the irrelevant operators are constructed using smeared links. While the use of such 
smearing allows for the use of highly improved definitions of the field strength tensor -F),^, we show 
that the standard 1-loop clover term with a mean field improved coefficient Csw is sufficient to 
remove the 0{a) errors, avoiding the need for non-perturbative tuning. This result enables efficient 
dynamical simulations in QCD with the FLIC fermion action 

PACS numbers: ll.15.Ha, 12.38.Gc 



I. INTRODUCTION 



The fat-link irrelevant clover (FLIC) fermion action [Tl 
is an efficient [2] Wilson-style nearest-neighbour lattice 
fermion action which incorporates both the thin gauge- 
field links of the Markov chain and fat links - links cre- 
ated via APE HYP or stout-Hnk d smear- 
ing. Through the use of fat links in the irrelevant oper- 
ators of the action, one achieves significant improvement 
in the chiral properties of the action reflected in a nar- 
rowing of the distribution of the critical Wilson mass Q . 
One also by-passes the fine-tuning problem typically en- 
countered in 0{a) improvement, as the use of fat links in 
both the irrelevant Wilson and clover terms suppresses 
the otherwise large renormalizations of the improvement 
coefficients. At the same time, short-distance physics is 
preserved completely in the action as the relevant oper- 
ators are constructed with thin links. 

Previous work established the good scaling proper- 
tics of the Fat Link Irrelevant Clover (FLIC) fermion ac- 
tion when a highly improved definition of the lattice field 
strength tensor F^n is used in the clover term. In this 
work we demonstrate that the use of the standard 1-loop 
definition of F^n with fat links in the clover term is suf- 
ficient to provide 0{a?) scaling for FLIC fermions. The 
1-loop variant has the advantage of maintaining a sim- 
ple force term when performing the molecular dynamics 
portion of a Hybrid Monte Carlo algorithm to generate 
dynamical configurations. 

In Sec. |TT] we highlight the essential features of the 
FLIC action with a particular emphasis on the various 
lattice field strength tensors used in the simulations. In 
addition the S'J7(3)-projection method used to create the 
fat links is outlined. In Sec. IIIII we describe the methods 
used to obtain an accurate scale determination on each 
lattice considered. Simulation parameters and scaling re- 
sults are presented in Sec. IIVI while correlation function 
properties are examined in Sec. [V] Conclusions are sum- 
marized in Sec. IVII 



II. FLIC FERMIONS 

The FLIC fermion action [1] is a variant of the clover 
action where the irrelevant operators are constructed us- 
ing smeared links 0, |3| , and mean field improvement [Tlj 
is performed. The key point is that short-distance physics 
is suppressed in the irrelevant operators. This allows an 
effective mean-field improved calculation of the clover co- 
efficient, required to match the Wilson and clover terms 
such that 0{a) errors are eliminated [To| . Further, the 
improved chiral properties of FLIC fermion action allow 
efficient access to the light quark regime [1] . 

The FLIC operator is given by 



(1) 



where the presence of fat (or smeared) links and/or mean 
field improvement has been indicated by the super- and 
subscripts. The mean field improved lattice gauge covari- 
ant derivative is defined by 

ymfi = XI T"^/^ {'Ut.ix) 5^+jx,y - Ulix - jl) 5^-f,^y) , 

. . (2) 

and likewise the (smeared link) lattice Laplacian is such 
that 

IXr. 

(3) 

We choose tr^i/ = f [7^,7iy]- For the clover term, one 
usually selects a standard one-loop F^^, 



F^A^)^--{C^Ax)-Cl{x)), 



-(-C/^,_^(a;) + C/_^,_^(a;)), 



(4) 



(5) 



where Ufj,^^{x) — Uf^{x) U„{x + fi) Uj^{x + v) Ul{x) is the 
elementary plaquette in the fx, v plane. However, with the 
use of fat links, one is also able to choose highly improved 
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definitions of Let C^^"(a;) correspond to the 

sum of the four m x n loops at the point x in the clover 
formation, and then define 



where mq and Uq are the mean links for the standard and 
fattened links. We calculate the mean link via the fourth 
root of the average plaquette 



prn X n 



{x) = 



Cfm X n 



(x)). 



(6) 



We can construct a 2-loop field strength tensor which is 
free of 0{a^) errors, 



1 , 



I ' 



or a 3-loop version which is free of O(a^) errors, 

p3L _ ^plxl ^p2x2 I J_p3x3 

^ iiv 2 20^" 90 '^'^ ' 



(7) 



(8) 



The smeared links in the FLIC action can be equally 
well constructed from standard APE smeared links, or 
the more novel stout link method [H. As the smeared 
links only appear in the irrelevant operators, the physics 
of the action are essentially independent of the choice of 
smearing method. The only requirement is that sufficient 
smearing is done such that the mean field improvement 
becomes an effective means of estimating the clover co- 
efficient Csw We typically find that four sweeps of APE 
smearing at a = 0.7 or four sweeps of stout smearing at 
p = 0.1 to be sufficient for lattices with a spacing between 
0.1 and 0.165 fm. 

In this work we use APE smeared links Ujl{x) con- 
structed from U^{x) by performing 4 smearing sweeps, 
where in each sweep we first perform an APE blocking 
step (at a — 0.7), 



V^^=\x)^{l-a)^ + 



(9) 



followed by a projection back into SU{3),uj:P (x) — 

'P(K'"'^(a;))- We follow the "unit-circle" projection 
method given in p^ . which allows for dynamical sim- 
ulations. The projection is defined by first performing a 
projection into U (3) 



U'{V) = V[V''V] 
followed by projection into SU (3) 

V{V) = , ^ 

^det U'{V) 



U'{V). 



(10) 



(11) 



It should be noted that the principal value of the cube 
root (being that with the largest real part) is the appro- 
priate branch of the cube root function to choose. As 
noted in (l3| this choice provides the mean link which is 
closest to unity. 

Mean field improvement is performed by making the 
replacements 



Uo 



K{x) 



(12) 



Wo = ( ^ ^eTrU^„{x) 



(13) 



a;, fi<u 



III. SCALE DETERMINATION 



The scale is determined using a 4-parameter ansatz 



V{r) = Vo + ar-e 



+ 1 



(14) 



as in Ref. [l4| . The tree- level lattice Coulomb term used 
in the ansatz is given by 











r 





2^ 



cos(k-r)i:»oo(0,k). (15) 



Here £'oo(0,k) comes from the tree- level gluon propa- 
gator for the appropriate gluon action. For the Wilson 
gluon action, we have at tree-level. 



^o"o'(0,k) 



4 ^ sm 



2 ^fJ- 



(16) 



where on a lattice with extents the allowed momenta 
are 



, ^ <n,,< 

' 2 ''-2 



(17) 



For the Liischer-Weisz gluon action, we have at tree- 
level. 



7^oo^(0,k)^4 5:fsin^^ + isin4^ 



(18) 



The lattice Coulomb term is constructed by calculating 
on large lattice volumes and then extrapolating to infinite 
volume. Explicitly, we choose L = 128 and L = 256 
and calculate [i] ^ for an spatial volume. On a finite 
volume, the Coulomb term takes the form [l5l | 



L-r 



1 
L 



O 



(19) 



In order to calculate the infinite volume tree-level lattice 
Coulomb term [i] , we extrapolate [i]^ linearly in to 



The tree-level lattice Coulomb term [i] for the Wilson 
and Liischer-Weisz gauge action is shown in Fig [T] The 
important finite lattice spacing artefacts are revealed at 
small r < 3a. The 0{a^) improvement in the Liischer- 
Weisz Coulomb term is also readily apparent. 
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FIG. 1: The (infinite volume) tree-level lattice Coulomb term 
for the Wilson and Liischer-Weisz (IMP) gauge action. 
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IV. SCALING RESULTS 

Calculations are performed on mean-field improved 
plaquette plus rectangle SU{3) Liischer-Weisz lattices. 
Lattice spacings determined using fits to Eq. ITil) above 
are given in Table [D 



/3 



H(fn 



a[ro](fm) 



4.80 
4.60 
4.53 
4.38 



0.096(1) 
0.120(1) 
0.132(1) 
0.164(1) 



0.088(1) 
0.113(1) 
0.124(1) 
0.152(1) 



TABLE I: The lattice spacing for pure Liischer-Weisz glue 
determined by the string tension ^/a = 440 MeV and the 
Sommer scale ro — 0.49 fm for various couplings (3. 

For each of the lattices we calculate quark propagators 
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Mntq 


Mpro 


FLIC-IL 


4.60 


2.278(26) 


3.347(33) 


2.638(30) 


3.875(39) 




4.53 


2.313(27) 


3.368(41) 


2.662(31) 


3.876(47) 




4.38 


2.299(21) 


3.323(32) 


2.688(25) 


3.886(38) 


FLIC-2L 


4.60 


2.347(26) 


3.394(33) 


2.717(30) 


3.929(39) 




4.53 


2.39(30) 


3.453(44) 


2.751(35) 


3.974(51) 




4.38 


2.41(24) 


3.450(35) 


2.818(28) 


4.034(41) 


FLIC-3L 


4.60 


2.365(30) 


3.461(37) 


2.738(34) 


4.006(43) 




4.53 


2.413(37) 


3.478(48) 


2.776(43) 


4.003(55) 




4.38 


2.435(27) 


3.474(38) 


2.847(32) 


4.062(44) 



TABLE II: Results for the A'^ and p masses on the three lat- 
tices, for the scale determined by the string tension a and the 
Sommer scale ro. 
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FIG. 2: The scaling of the and p masses for various quark 
actions in the quenched approximation according to the string 
tension (upper) and the Sommer scale (lower). 



using the FLIC fermion action with a 1, 2 and 3-loop 
clover term as described in Sec.|lll The tt, p and N masses 
are then calculated and interpolated to a 7r/p mass ratio 
of 0.7, shown in Table HH 

Scaling results are presented in Fig. [51 The lines of 
fit are extrapolations in constrained to pass through 
the single point at the continuum limit. The lines for 
non-perturbatively improved clover and all FLIC actions 
are straight, indicating O(a^) scaling, that is the effective 
elimination of 0{a) errors. 
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Thus, 1-, 2- and 3-loop fat- link formulations of i^^jy in 
the FLIC fermion action all provide 0{a) improvement 
as expected. The different formulations differ at the level 
of O(a^). Remarkably, the 1-loop action is actually the 
preferred action. Firstly, it is the cheapest to perform 
molecular dynamics with, which is important for Hybrid 
Monte Carlo dynamical simulations. Secondly it has the 
smallest residual 0{a^) errors in the quantities we have 
studied here. We'll also see that correlation functions 
have smaller fluctuations. 



V. CORRELATION FUNCTIONS 

Finally, we compare the p— meson correlation function 
on the fine (3 = 4.60 and coarse /? = 4.53 lattices at ap- 
proximately matched pion masses for the three different 
FLIC actions. The source is at time slice 8. 

The effective mass plots are given in Figure [31 The 
main effect that we observe is that as the Euclidean time 
index progresses into the latter half of the lattice, the 
1-loop FLIC correlators show reduced fluctuations and 
reduced error bars when compared with the 2-loop and 
3-loop FLIC results. The difference is particularly ob- 
servable on the coarser (3 = 4.53 lattice. We understand 
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this to be due to the 1-loop action having a more lo- 
cal field strength f^^(x) than the 2- and 3-loop actions 
making it less susceptible to large fluctuations. 

VI. CONCLUSIONS 

We have examined the role of improvement in the lat- 
tice field strength tensor of the FLIC fermion action, Our 
results demonstrate that the standard 1-loop choice of for 
the lattice clover term in the FLIC fermion action pro- 
vides 0{a?) scahng. 

Remarkably the 1-loop action provides results that 
are preferable to those obtained from the 2-loop 0{a?)- 
improved lattice field strength tensor or those obtained 
from the the 3-loop C'(a'*)-improved definition. The 1- 
loop results provide 

1. Smaller residual 0{a^) errors, 

2. Stable hadron correlators with reduced fluctua- 



tions, 

3. Smaller statistical uncertainties, and 

4. A more efficient action suitable for dynamical 
fermion simulations. 

This result enables efficient and effective dynamical 
QCD simulations with FLIC fermions. Simulations are 
currently under way. 
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